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Abstract 

We compute in order Ug the nonrelativistic QCD (NRQCD) sliort-distance coefficients that match 
quark-antiquark operators of all orders in the heavy-quark velocity v to the electromagnetic current. 

^; 

^ . We employ a new method to compute the one-loop NRQCD contribution to the matching condition. 

The new method uses full-QCD expressions as a starting point to obtain the NRQCD contribution, 

^ ' thus greatly streamlining the calculation. Our results show that, under a mild constraint on 

cn . the NRQCD operator matrix elements, the NRQCD velocity expansion for the quark-antiquark- 

\0 ■ 

^^ \ operator contributions to the electromagnetic current converges. The velocity expansion converges 

l> 

^^ ' rapidly for approximate J/'0 operator matrix elements. 
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I. INTRODUCTION 

The electromagnetic decays of quarkonia through a single virtual photon have played 
an important role in the experimental and theoretical development of quarkonium physics. 
On the experimental side, the decays of charge-conjugation-odd quarkonium states to a 
lepton pair provide unique signals for the detection of those states. On the theoretical 
side, the decays of ^5*1 quarkonium states to a lepton pair allow one to determine one of 
the fundamental parameters of the heavy-quark-antiquark (QQ) bound state, namely, the 
square of the wave function at the origin. (See, for example, Ref. [1].) The square of the wave 
function at the origin enters into many calculations of quarkonium decay and production 
rates. 

The expression for the ^5*1 quarkonium decay rate into a lepton pair at leading order in 
the QCD coupling ag and at leading order in v, the Q or Q velocity in the quarkonium 
rest frame, has been known since the first discovery of quarkonium and is based on the 
Van Royen-Weisskopf formula [2] of quantum electrodynamics. The order-ttg corrections to 
this formula at leading order in v were calculated in Refs. [3, 4]. The order-a^ relativistic 
corrections at relative orders f^ and v^ were calculated in Refs. [5, 6], respectively. Order- 
ctg corrections to the decay rate were calculated in Refs. [7, 8]. The correction to the 
electromagnetic current of a quarkonium at relative order agv'^ was calculated in Ref. [9]. 

In this paper, we calculate relativistic corrections to the quarkonium electromagnetic 
current at order a.^- We carry out our calculation in the context of nonrelativistic QCD 
(NRQCD) [5]. We obtain closed- form expressions whose Taylor-series expansions in v give 
the short-distance coefficients for the NRQCD QQ operators, of all orders in v, that match 
to the electromagnetic current. We do not consider QQ operators that contain gauge fields. 
Therefore, our operators are not gauge invariant, and we evaluate their matrix elements 
in the Coulomb gauge. In the Coulomb gauge, QQ operators involving gauge fields first 
contribute at relative order f*^. Our results confirm the calculation at relative order agv'^ in 
Ref. [9]. Since the corrections at relative order agv"^ are not very significant at the current 
level of precision of calculations of '^S'l quarkonium electromagnetic decay rates, we do not 
expect the order-a^ corrections at still higher orders in v to be important numerically. 

We present our calculation primarily as a demonstration of a new method for comput- 
ing the one-loop NRQCD contribution that enters into the matching of NRQCD to full 



QCD. The direct computation of one-loop NRQCD expressions to all orders in v would 
be a formidable task, in that it would require knowledge of the NRQCD interactions and 
electromagnetic-current operators, their Born-level short-distance coefficients, and their 
Feynman rules to all orders in v. Instead of following the direct NRQCD approach, we 
note that NRQCD through infinite order in v is equivalent to QCD, but with the inter- 
actions rearranged in an expansion in powers of v. Therefore, we can obtain the one-loop 
NRQCD contribution by starting from full-QCD expressions and expanding integrands in 
powers of momenta divided by the heavy-quark mass m before we carry out the dimen- 
sional regularization. In dimensional regularization, this method is related to the method 
of regions [10]. We explain this relationship in Sec. II. The method of regions has been 
used previously at leading order in v to compute NRQCD short-distance coefficients from 
full-QCD expressions. (See for example, Ref. [8].) 

Our results show that, under a mild constraint on the NRQCD operator matrix elements, 
the NRQCD velocity expansion for the quark-antiquark-operator contributions to the elec- 
tromagnetic current converges. The velocity expansion converges rapidly for approximate 
J/ip operator matrix elements. 

The remainder of this paper is organized as follows. In Sec. II we discuss the one-loop 
matching of NRQCD to QCD at all orders in v. We define the notation that we use to 
describe the kinematics of the calculation in Sec. III. Section IV contains detailed formulas 
for the NRQCD QQ short-distance coefficients. In Sec. V we compute the one-loop QCD 
corrections to the electromagnetic current, while in Sec. VI we use our new method to 
compute the one-loop NRQCD corrections to the electromagnetic current. We give analytic 
and numerical results for the short-distance coefficients in Sec. VII, present a formula that 
resums a class of relativistic corrections to all orders in f , and discuss the convergence of 
the velocity expansion. Our conclusions are given in Sec. VIII. The Appendices contain 
compilations of integrals and identities that are useful in the calculation. 

II. MATCHING TO ALL ORDERS IN v 

We define the hadronic part of the quarkonium electromagnetic decay amplitude A^^ as 

{-ieeQ)iA%={Q\JUH), (1) 



where H is the quarkonium, e is the electromagnetic charge, eg is the heavy-quark charge, 
and Jgj^ is the heavy-quark electromagnetic current: 

J^M = (-^eeQ)^7^^. (2) 

Here, ip is the heavy-quark Dirac field, and 7^^ is a Dirac matrix. 

In the quarkonium rest frame, iA^j^ = because of conservation of the electromagnetic 
current. According to NRQCD factorization [5], we can write the spatial components iA^jj 
as 

zA'h = v/2^5Zc„(0|O:ji/), (3) 

n 

where the c„ are short-distance coefficients, the (9^ are NRQCD operators, and itlh is the 
quarkonium mass. We regulate the operator matrix elements in Eq. (3) dimensionally in 



d = A — 2e dimensions. The factor ^Imu on the right side of Eq. (3) appears because the 
NRQCD operator matrix elements have nonrelativistic normalization, while we choose the 
amplitude on the left side of Eq. (3) to have relativistic normalization for the quarkonium 
E. 

The aim of this paper is to calculate the short-distance coefficients c„ that correspond 
to QQ color-singlet operators in order a\. We can determine these c„ by making use of 
a matching equation that is the statement of NRQCD factorization for perturbative QQ 
color-singlet states: 

^^bOi = E^"<0|^"l^'5i), (4) 

n 

where the subscript 1 indicates a color-singlet state. Throughout this paper, we suppress the 
factor ^JWc that comes from the implicit color trace in iA^^ , where A*";, = 3 is the number 
of colors. Through order a], the matching equation is 

tX^l+t^^l = $:(c(°) +c«)(0|O;|QQi)(°) + 5^cW(0|O^|Qgi)«, (5) 

n n 

where the superscripts (0) and (1) indicate the order in a^. In the first sum in Eq. (5), 
only color-singlet QQ operators contribute, while in the last sum, additional operators can 
contribute if they mix into color-singlet QQ operators under one-loop corrections. 
We define the quantity 



which is the expansion of i-^nQ in powers of g/m, where q is half the relative momentum 
of the heavy quark and heavy antiquark. At order q;°, the matching equation (5) yields 

^<i = E^«°^(0|O;|QQ,)(°), (7) 

n 

from which the cii can be determined. The Cn have been computed previously in Ref . [1] . 
At order aj, the matching equation (5) yields 

<A = E -i^^ (^\On\QQ.r' + KA] ^^^^^ , (8) 



iA 



i(l) 



NRQCD 



from which the c„ can be computed. We compute the quantities iA^J^ and 
in Sees. V and VI, respectively. 
The quantity 

n 

would be formidable to calculate directly in NRQCD because it involves operators and 
interactions of all orders in v. Rather than carry out such a direct calculation, we take a 
new approach. We note that, by construction, NRQCD reproduces all of the interactions in 
full QCD, but with those interactions reorganized in an expansion in powers of v. Therefore, 



we can obtain 



iA 



from the expression for iA^^I by expanding the integrand in 

. NRQCD '^'^1 

powers of the momentum divided by m. Before making this expansion, we carry out the 
integration over the temporal component of the loop momentum, using contour integration. 
This procedure establishes the scale of the temporal component of the loop momentum, 
which varies from contribution to contribution, and it avoids the generation of ill-defined 
pinch singularities that can arise when one expands the Q and Q propagators prematurely 
in powers of the momentum. We expand the integrand in powers of both the external 
momenta divided by m and the loop momenta divided by m. We then regulate the integrals 
dimensionally, setting scaleless, power-divergent integrals equal to zero. Ultimately, we 
renormalize ultraviolet divergences according to the MS prescription. 

The procedure of expanding both external and loop momenta in power series before 
regulating dimensionally was first utilized in the Appendix of Ref. [5]. The rationale for it 
was discussed in Refs. [9, 11]. This procedure amounts to the prescription that infrared- 
finite contributions that arise from loop momenta in the vicinity of zero are kept in the 



short- distance coefficients [11].^ 

If one uses dimensional regularization, then the quantity ^n'^^i \0\Ol^\QQi)^^'^ corre- 
sponds to the contribution from the hard region in the method of regions [10], while the 
quantity 



^•^QQi 



corresponds to the sum of the contributions from the potential, soft, 

NRQCD 

and ultrasoft regions, i.e., the contribution from the small-loop-momentum region. In the 



method of regions, it is assumed that there are no contributions from the region in which 
the temporal component of the gluon momentum is of order m, but the spatial component 
of the gluon momentum is of order mv. As we shall see explicitly in our calculation, this 
assumption is justified because the contribution from this region of integration vanishes in di- 
mensional regularization. We note, however, that the contribution from this region does not 
vanish in the case of a hard-cutoff regulator. One potentially useful feature of the approach 
that we present here is that it can be applied in the case of a hard cutoff, such as lattice 
regularization, while the method of regions is applicable only in dimensional regularization. 
In the method of regions, one can compute the contribution from the hard region directly, 
rather than computing it, as we do, by subtracting the small-loop-momentum contribution 
from the full-QCD contribution. As we shall explain later, there may be advantages to our 
indirect procedure in calculating the hard contribution to all orders in v. 

III. KINEMATICS 

Before proceeding to write explicit formulas for the short-distance coefficients, let us 
define some notation for the kinematics of the heavy-quark electromagnetic vertex. We 
take pi and p2 to be the momenta of the incoming heavy quark Q and heavy antiquark Q, 
respectively, pi and p2 can be expressed as linear combinations of their average p and half 
their difference q: 

Pi = P + q, (10a) 

P2 = p- q- (10b) 



^ In the case of hard-cutofF regularization, such as lattice regularization, the expansion in powers of loop mo- 
mentum divided by m would be uniformly convergent and would yield the same result as the unexpanded 
expression. 



In the QQ rest frame, the momenta are given by 

Pi = {E,q), (11a) 

P2 = {E,-q), (lib) 

p=(E,0), (lie) 

g=(0,q), (lid) 

where E = ^/m? + q^. The quark Q and antiquark Q are on their mass shells: Pi = P2 = ''^^• 
For later use, it is convenient to define a parameter 

^-|. (12) 



which is related to the velocity 



We can write S in terms of v as 



VTTv 

E"^ and q^ are expressed in terms of m and 6 as 

,2 



t; = M. (13) 

m 



(14) 



n2 m 



9 = 137^- (15b) 



IV. FORMULAS FOR THE SHORT-DISTANCE COEFFICIENTS 

Now let us make use of the matching conditions (7) and (8) to compute the short- distance 
coefficients for the specific color-singlet QQ operators that we consider in this paper. These 
operators are 

0\n = xH-fV)2VV, (16a) 

O'sn = X^-f V)2"-2(-fV0(-f V) ■ CT^, (16b) 

where ip is the Pauli spinor field that annihilates a heavy quark, x^ is the Pauli spinor 
field that annihilates a heavy antiquark, and o"* is a Pauli matrix. Our operators contain 
ordinary derivatives, rather than covariant derivatives. Therefore, our operators are not 
gauge invariant, and we evaluate their matrix elements in the Coulomb gauge. We do not 



consider QQ operators involving the gauge fields, which first contribute at relative order v'^. 
We note that (9^„ can be decomposed into a linear combination of the S'-wave operator 0\^ 
and the D-wave operator C^„: 

1 



-S" ~ fl — 1 ^" ~'~ Dni 



(17) 



where C})„ is defined by 



o'on = x\-wr-' 



;-m(-iv)-a 



1 



-yva 



2_i 



i,. 



d-V 2 

In the basis of operators 0\^ and 0^„, the matching conditions (7) and (8) become 



(18) 



^<A = E«"^(0|^VIQOl)^°^ +E^n^(0|^BnlWl) 



(0) 



^<A = E«"^(0|^VIQQl)^°^ +E^n^(0|^BnlWl)^°^ + 



M' 



;(i) 



NRQCD 



(19a) 
(19b) 



where a„ and 6„ are the corresponding short-distance coefficients. A similar equation holds 
in the basis 0\^ and 0^„, where the associated short-distance coefficients are 



a„ + 



Jji '-^n 



'^n "n; 



d-l 



(20a) 
(20b) 



respectively.^ 

The QQ matrix elements in Eq. (19) are 

(o|oviWi)^°) 
(o|oi,jggi)(°) 






(21a) 
(21b) 



If we replace the ordinary derivatives V with covariant derivatives D in an S'-wave operator 0\^, then we 
obtain one of the conventional gauge-invariant S'-wave NRQCD operators. Because the squared covariant 
derivatives {D)"^ commute with themselves, the substitution of covariant derivatives for ordinary deriva- 
tives leads to a unique S'-wave operator at each order n. Therefore, the S-wave short-distance coefficients 
s„ that we compute are also the short-distance coefficients of the S-wave operator in which ordinary 
derivatives have been replaced with covariant derivatives. In the case of the D-wave operators O^j^, the 
replacement of ordinary derivatives with covariant derivatives does not lead to a unique operator because 
{Dy and {Dy do not commute. Therefore, each of the D-wave short-distance coefficients dn that we 
compute is the sum of the short-distance coefficients for the various operators at order n that can be 
constructed from covariant derivatives. 



where ^ and rj are two-component spinors. In order to maintain consistency with our calcu- 
lations in full QCD, we have taken the QQ states to have nonrelativistic normalization and 
we have suppressed the factor y/W^ that comes from the color trace. 
Because of current conservation, the most general form of iA\/^ is 

' ^ QQi 

tA^QQ^=v{p2){G^' + Hq^)u{p,), (22) 

where 

G = Zq{1 + A). (23) 

Zq is the fermion wave-function renormalization, and A is the multiplicative correction to 
the fermion electromagnetic vertex. Similarly, 



A' 



NRQCD 



^(P2)(G'nrqcd7* + H^mcDq')u{pi), 



where 



G 



[Zq 



NRQCD — [^QJ NRQCD 



:i + A 



NRQCD J 



Using nonrelativistic normalization for the spinors u and v, we obtain 



v{p2)Yu{pi) = V 0-'^ 






q'-v{p2)u{pi 



q'i]^q ■ (Tj 

E 



Then, 



Similarly, 



iA'qq^ = Gr^^a^i 



G 



H 



_E{E + m) E_ 



q'ri^q ■ a^ 



A, 



QQi 



NRQCD 



Gnrqcd^V*^ 



t^V 



G 



NRQCD 



+ 



H, 



NRQCD 



_E{E + m) E 



gVq ■ 0-6 



(24) 
(25) 

(26a) 
(26b) 

(27) 
(28) 



Using the matching condition (19a) and Eqs. (21) and (27), we obtain the short- distance 
coefficients at order a2'- 



al°) = 1 ^ ^ 



n\ \dq 



cw 



5, 



q2=0 



hf = df 



d 



(n-1)! ydq"^ 

n-l 



nO) 



n-l 



(29a) 



G-Co) /j(o) 



_E{E + m) E 



q2=0 



d 



(n — 1)! \dq^ 



_E{E + m)_ 



q2=0 



n n 



(29b) 
(29c) 



Using the matching condition (19b) and Eqs. (21), (27), and (28), we obtain the short 
distance coefficients at order a]: 

I f d 



a(') = - 
n! 






1 



(n 



.(1) 



ai^) + 



1 



d 



-1) 
6«, 



_d_ 
9^ 



71-1 



E{E + m) 



+ 



A/J(i) 



E 



where 



AAUr — Ur "-^NRQCD' 

AiJ(i) - iJ(i) - iJ^^^ 



(30a) 
(30b) 
(30c) 

(31a) 
(31b) 



The infrared divergences in Gj^jj^q^^j^ and -f^NROCD cancel in AG*^^^ and AH'^^^ because 
NRQCD reproduces full QCD in the infrared region. The one-loop NRQCD matrix ele- 
ments in Gp^f{^Q,2:D contain ultraviolet divergences, which we renormalize according to the 
MS prescription. The quantity -f^NRQCD i^ ^^^ of ultraviolet divergences. The quantities 
A and Zq also contain ultraviolet divergences. However, because of the usual cancellation 
between the vertex and fermion-wave-function renormalizations, G^^^ is free of ultraviolet 
divergences. H^-^^ is also free of ultraviolet divergences. Carrying out the renormalization, 
we have 



r « J MS ^ ~ 



MS 



n J MS 



d 



1)! \dq 



n-l 



MS _|_ ^-" 



E{E + m) 



E 



L"'" J MS In J MS ^ 3 In J MS 



(32a) 
(32b) 
(32c) 



In deriving the expression for 



,(1) 



MS 



we have used the fact that, in minimal subtraction. 



one removes the 1/e pole times the order-a^ rf-dimensional matrix element. Hence, a term 
proportional to {d — l)^^e^^ is subtracted in Eq. (30c) in carrying out the renormalization. 



V. QCD CORRECTIONS 



In this section, we calculate the QCD corrections to the heavy-quark electromagnetic 



current. That is, we compute iA 



i(i) 



10 



A. Vertex Correction 



In the Feynman gauge, the vertex correction to the electromagnetic current is given by 



AM = -rg^^Cp 



2r< I ^iP2hai-h + 1^ + "^)7''(A + ^ + m)-f"'u{pi) 



Ik D0D1D2 

where g"^ = Anas is the strong couphng, Cp = {N^ — l)/{2Nc) = 4/3, and 



(33) 



/i 



2e 



Do = k"^ + ie, 

Di = P + 2k-pi +ie, 

D2 = P -2k-p2 + ie. 



(34a) 
(34b) 
(34c) 
(34d) 



/i is the renormahzation scale. The loop momentum k is chosen to be the gluon momentum. 
By making use of Eq. (10) and applying the equations of motion. 



v{p2)i)u{pi) = 0, 
v{p2)iu{pi) = mv{p2)u{pi), 



(35a) 
(35b) 



we find that Eq. (33) can be written as 

4 DoDrD2 



A^ = -iglCp I ^ ^ ^ 'j\:p2j 



■v{p2 



{d - 2)A;2 - 4(2p2 -m^)+^k-q 



r 



+Amk^ - 8g^^ + 2(2 - rf)P^ Wpi). (36) 



Tensor reductions of the integrals in Eq. (36) are given in Appendix A. The result is 



AM = -rglCFv{p2) 



{d - 2) Ji - A{2p' - m')J2 + 4J3 + 2(2 - rf) J4 



Y 



2mp^ 2mq^ 



-J' 



p. q 

where the integrals Ji are defined by 

J,= 



J3 + 2(2 - d)m \ Jf 



Ni 



P^ J 

pZgZ 



Mpi). (37) 



DoDiD2 



(38) 



11 



and 






1, 

2k -q, 

1 
d-2 
2k -p, 

1 
d-2 
k ■ pk ■ q 



e- 



{k-pf {k-qf 



pz 



-e + 



{k-pf 

p2 



+ {d-l 



{k-q) 



21 



(39a) 
(39b) 
(39c) 
(39d) 
(39e) 
(39f) 
(39g) 



The integrals J1-J7 are evaluated in Appendix B. The results are tabulated in Eq. (B12). 
We note that J5 and Jj vanish, as is required by conservation of electromagnetic current in 
Eq. (37). 

Writing the vertex correction as A'^ = v{p2){h^'^^ + Hq'^)u{pi), we have 



A 



-igjCt 



{d - 2) Ji - 4(2/ - m^) J2 + 4J3 + 2(2 - rf) J4 



^ + log ^ , + 2 1 + 5^)L{5) — + log ^ ^ 

4vr I euv ^ V^iR "^ 



iJ = -tglC, 



A{l + 5^)K{5) + {l + 5^) 



2m 2(2 -d)m ' 
^Js H 5 Jq 



- 66^L{6) 
352 



Tc in f 1 TCfi e ^ 



asCp 1 - 5^ 



2L{6) - 



ZTT 



(40a) 
(40b) 



Air m 

where the subscripts on 1/e denote the origins of the divergences and 7^, is the Euler- 
Mascheroni constant. The functions L{6) and K{6) are given by 



^(^)-^'-(l^ 



^^■<*) = i 



Sp 



26 
TT6 



Sp 



26 



1-6 



(41a) 
(41b) 



where Sp is the Spence function: 



Sp(a;) 



log(l-t) 



dt. 



(42) 



In Eq. (40), we have neglected terms of order e^ and higher. In the remainder of this paper, 
we drop such higher-order terms. 



12 



B. Wave-function Renormalization 



The heavy-quark wave- function renormahzation Zq, evaluated in dimensional regulariza 
tion, is given in Ref. [12]: 

agCp 



Zq = 1 + 



An 



1 2 47r/x^e-^ 

3 log 4 

euv eiR m^ 



(43) 



C. Summary of QCD results 



By making use of Eqs. (23), (40), and (43), we find that G and H are given by 



G = 1 



asGp 
An 



(2 [(1 + 5^)1(5) -l](— + log ^""^^ ^ ] + 66^1(6) - 4(1 + 6^)Ki5) 



H 



-4 + (1 + 5^ 
a^Gp 1 - 52 r 



''tt^ in / 1 



, log 1 

S S VeiR q' 1 + 52 



2iW - J 



An m 
Expanding Eq. (27) through order v"^, using Eq. (44), we obtain 



^^QO. = ^'^'^ 



asGp (89/1 
l + ^-^l-v^ ' 



47r 1^3 VeiR 
,. 3t;2 



log 



An 



2m? 



m^ 

vr^ in f I 
V V VeiR 

TT^ ivr / 1 



log 



9 



t' t; VeiR 



log 



TT/i^e 'fe 



(44a) 
(44b) 



?>inv 



(45) 



Equation (45) agrees with Eq. (4.16) of Ref. 



VI. NRQCD CORRECTIONS 



In this section, we calculate the NRQCD corrections to the heavy-quark electromagnetic 



current. That is, we compute 



iA 



i{i) 



. In order to demonstrate our method for 



QQiJnrqcd 
calculating these corrections from fuU-QCD expressions, we present the calculation in some 

detail. 



13 



Divergent integrals are regulated using dimensional regularization, with d = 4 — 2e. We 
define the following notation for the loop integrals in rf — 1 dimensions: 



We also define jj}^ and jj}^, which have the same meaning as J^ and J^^, except that it is 
understood for jr^ that one carries out the k^ integration first, and it is understood for both 
jj} and j/ that one expands the integrand in powers of the momenta divided by m. 



A. Vertex Correction 



Now, we calculate the NRQCD vertex correction to the electromagnetic current. From 
Eq. (36), we see that the vertex correction is given by 



A' 



^^Q^^ = -^^'^^:^^^"(^^ 



{d - 2)A;2 - 4(2/ - m^) + 8k ■ q 



Y 



+4mk' - 8q% + 2(2 - d)¥J/i \u{pi). (47) 



The vertex correction (47) can be written as 



Anrqcd = -igiCFv{p2) (rf - 2)Si - A(2p' - m')S2 + Sq.Sl^ 



r 



+2 



2mSl - A-i^St^q' + (2 - rf)7^5r U(pi 



where 






Si; 



D,D.' 



*J4 



fkDo 

r. k^'k" 



(48) 

(49a) 
(49b) 
(49c) 
(49d) 



The factors in the denominator of the integrands are defined in Eq. (34). In the QQ rest 
frame, the factors Di in Eq. (34) are 



Do = {ky -k' + ie = (r - |fc| + ie)(r + |fc| - ie), 

D^ = (A;° + Ey-A^ + ie = {k^ + A + E-ie){k^-A + E + ie) 

D2 = (A;° - Ef - A^ + ie = {k^ + A - E - ie){k° - A - E + ie) 



(50a) 
(50b) 
(50c) 



14 



where A is defined by 

A = ^m^ + (fc + qf. (51) 

The following are identities that we use frequently: 

A + E ^ ' 

t^^ -{E±\k\f = ^2\k\{E^q-k), (52b) 

where a = a/\a\ for any spatial vector a. We first evaluate the k^ integral by contour 
integration, closing the contour in the upper half-plane in every case. The contributions 
from the poles in the gluon, quark, and antiquark propagators are defined as Sig, SiQ, and 
SiQ, respectively. Certain integrals that we use frequently are tabulated in Appendix C. 

We note that the contributions S^q correspond to the potential region in the method of 
regions, and the contributions Sig correspond to the soft and ultrasoft regions in the method 
of regions [10]. The contributions SiQ correspond to a region of integration in which the 
temporal component of the gluon momentum is of order m, but the spatial component of 
the gluon momentum is of order mv. As we have mentioned, in the method of regions it 
is assumed that this region of integration does not contribute [10]. We shall see explicitly 
in the calculations that follow that this assumption is justified because the contributions 
from this region of integration consist of scaleless, power-divergent integrals, which vanish 
in dimensional regularization. In the case of a hard-cutoff regulator these contributions do 
not vanish, and they must be included in the calculation of the NRQCD corrections. 

1. Si 

The integral 5*1 is the sum of two contributions: 5*1 = Siq + Siq. 

By making use of Eq. (50), we evaluate the k^ integral. The contribution from the quark 
pole is 



^"=""8l' 



Expanding 1/A and 1/(A + E) in Eq. (53) in powers of {k + qY/m'^, we find that all of the 
terms in the expansion are scaleless, power-divergent integrals. Hence, 

SiQ = 0. (54) 
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The contribution from the antiquark pole is 
We use the identity (52a) to reduce the integrand in Eq. (55) to the following form: 



c __ ^ 



fX'^DwT^T^e- (^^) 



Expanding 1/A in Eq. (56) in powers of (A; + qf' jvn? ^ we find that the expansion brings in 
additional factors of (fc + q)^. In each additional factor, only the term q^ survives, as the 
terms k^ + 2k ■ q lead to scaleless, power-divergent integrals, which vanish. As a result, we 
can replace A with E in Eq. (56). Hence, Siq is proportional to an elementary integral ni, 
which is defined in Eq. (C3a): 

Using Eqs. (12), (54), and (57), we obtain 

5, = ^ ^.S. (58) 

2. So 



The integral S2 is the sum of three contributions: S2 = S2g + S2Q + S2Q. 

By making use of Eq. (50), we evaluate the k^ integral. The gluon-pole contribution is 

1 



s. 



2 Jfc |fc|[A2 -{E+ \k\)^ - ie][A^ -{E- \k\)^ - ie] 
' ' ' (69) 
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8 A|fc|3[E2-(q.fc)2] 

where we have used the identity (52b). Making use of Eq. (C4b), we find that S2g is 
proportional to no in Eq. (C2): 

Using Eqs. (15) and (C2), we express S2g in terms of m and 6 as 

327r2m2 Veuv eiR/ 2o \1 — I 
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The contribution from the quark pole is 






A{A + E)[{A + Ey - k^ + ie] 



n=0 "^'^ 



A{A + E) 



2n+3' 



(62) 



Now we expand 1/A and 1/(A + E) in Eq. (62) in powers of {k + qY /irP'. All of the terms 
in the expansion yield scaleless, power-divergent integrals, which vanish. Therefore, we have 



S- 



2Q 



0. 



The contribution from the antiquark pole is 



S. 



2Q 



'se£] 



1 



A(A -E- ie)[k'^ - (A - E)^ - ie] ' 
If we use the relation (52a), we obtain 

E\ 1 



S. 



2Q 



-M 



1 + 



A 



fc2(fc2 + 2k-q-i€) 



1 1_ / fc2+2fcg 

^ fc2 \ A+E 



(63) 



(64) 



(65) 



The denominator factor in the brackets can be expanded to give 

1 ^ (fc2 + 2fc ■ g)2"-i 



S. 



2Q 



MH 



\2n 



(66) 



k\k^ + 2k ■ q - is) ^ fc2n+2(^ + ^)2 

Now we expand E/A and 1/(A + E) in powers of {k + qY /m^. The expansion brings in 
additional factors of {k + q)^ in each term in the integrand of Eq. (66). In each additional 
factor (fc + g)^, only the term q^ survives, as the terms k^ + 2k ■ q lead to scaleless, power- 
divergent integrals. Therefore, we can replace A in Eq. (66) with E. Furthermore, in the 
numerator of the second term in brackets in Eq. (66), only the term {2k ■ qr)^"^^ survives, 
as the other terms lead to scaleless, power-divergent integrals. Then, we have 



s. 



2Q 



M 



^i^(fc.q) 



(67) 



k'^ik^ + 2k ■ q - ie) 2^^ fc2n+2^2n 

^ ' n=l 

The term proportional to {k ■ qr)^""^ yields a scaleless, logarithmically divergent integral. 
However, this integral vanishes because the integrand is an odd function of k. Thus, only 
the first term in the brackets in Eq. (67) survives, and we find that 



S20 = -TTT^2 



1 



4E 



647[ E\q\ \eiR 



1 vr/i^e '4: 

— + log — - — 



tTT 



(68) 
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where n2 is defined in Eq. (C3b). 

Making use of Eqs. (15), (61), (63), and (68), we obtain 

i 1-6 



S,. 



(47r)2 4m2 



2L{6){ — - — )-- + -(— + log ^ 



euv eiR/ S 5 Vcir q^ 



(69) 



wfiere L{6) is defined in Eq. (41a). 



3. S^ 

Tlie integral 5*3 is the sum of three contributions: S^ = Sl^g + S^q + S^^. 

We first evaluate S^. The integral of S^ over A;° is identical to the integral of 5*2 over 
k^ except that, in S^, the result contains an additional factor of k^ evaluated at the gluon, 
quark, or antiquark pole. Thus, by making use of Eqs. (52a), (59), (62), and (66), we obtain 



1 



SL = — / ^ — , (70a) 

5*3 is a scaleless, power-divergent integral, which vanishes. In S^q and S^^ we expand 1/A 
and 1/(A + E) in powers of {k + qY /rn^. We find that every term in the expansions leads 
to a scaleless, power-divergent integral, which vanishes. Hence, 

St = 0. (71) 

Next we compute the spatial components S^. The integral of S^ over k^ is identical to 
the integral of 5*2 over k^ except that, in S^, the result contains an additional factor of /c*. 
Thus, by making use of Eqs. (59) and (62), we obtain 

SL = - [ ^, (72a) 

' 8jk\k\^[E^-{q-ky] 



n=0 ''^ 



iU'in 



sh = -T^y.4r .,/''^„,.. - (T2b) 



5*3 is a scaleless, power-divergent integral, which vanishes. Expanding 1/A and 1/(A + E) 
in 5*3^ in powers of (fc + q)^/m^, we also obtain only scaleless, power-divergent integrals, 
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which vanish. If we multiply the second term in brackets in Eq. (67) by k^, we obtain only 
scaleless, power-divergent integrals. Hence, 

if k^ 

After making a standard reduction of the tensor integral in Eq. (73) to a scalar integral, we 

obtain 

i ^ iq' r (fc^ + 2fc ■ q) - fc^ ^ iq' ^ 1 g^ 

^^ 8Eq^ J^k^{k^ + 2k-q-ie) SEq^""^ 32n E\q\' ^ ' 

where rii is defined in Eq. (C3a) and we have discarded scaleless, power-divergent integrals. 



Hence, 



1 0* 
SI = — —4r-r (75) 

^ 327rE|g| ^ ' 



Writing our results in Eqs. (71) and (75) in covariant form, we obtain 



^3 " (4vr)2 2m2 5 ^ ' ^^^> 

where we have made use of Eq. (15) to express E and \q\ in terms of 5. 

The integral Sj^'^ is the sum of three contributions: S^'^ = S^^ + S^q + S^q. 

We first evaluate 6*4°. The integral of Sf^ over k^ is identical to the integral of S^ over 
k^ except that, in 6*4°, the result contains an additional factor of k^ evaluated at the gluon, 
quark, or antiquark pole. Thus, by making use of Eqs. (52a) and (70), we obtain 

1 



'" SJk\k\[E^-{q-ky] ^ ' 



4Q 

n=0 

°" " (fc2 + 2fc ■ q)2«+l 



qOO 



i ^ r (k^ _^ 2k ■ 0)^"+^ 

"SE ^ ;5fc fc2n+2^(^ + ^)2n+l " (77c) 

Every integral in Eq. (77) is a scaleless, power-divergent integral. Hence, 

Sf^ = 0. (78) 

Next we compute 5*4*. The integral of 5*4* over k^ is identical to the integral of 5*3 over 
A;° except that, in S^\ the result contains an additional factor of k\ Thus, by making use 
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of Eq. (70), we obtain 



COJ 






cOi_ 



t 



SJkk^[E^-{q-ky 



^ r. k\k^ + 2k-qf'' 



I 

SE 



fc2n+2^(^ + E) 



(79a) 
(79b) 
(79c) 



5*4* is a scaleless, power- divergent integral, which vanishes. S^q and 5*°)^ also vanish, once 
we expand 1/A and 1/(A + E) in powers of {k + qY /m^. Thus, 



^f 



0. 



(80) 



Finally, we evaluate the integrals S]l . The integral of S]^ over k^ is identical to the integral 
of 5*3 over A;° except that, in S'4"', the result contains an additional factor of y . Thus, by 
making use of Eqs. (72) and (73), we obtain 



nij _ ^ / 2_J 

^' ~ ^Jk\k\^[E^-{q-kf]' 



qij _ _ 

^« - 8E 



qij 
•^4(9 



00 „ 

n=0 •^'^ 



k'^k 



2n 



A(A + E)2'^+3' 

k'y 



^la) 



(81b) 



fSlc) 



4E J^k'^{k'^ + 2k-q-ie)' 

S'Z is a scaleless, power-divergent integral, which vanishes. 5*4^ also vanishes, once we 
expand 1/A and 1/(A + E) in powers of {k + qY /m^ . The tensor integral S^^q in Eq. (81c) 
must be a linear combination of the two symmetric tensors 5"^^ and q^q^ . By contracting 
these tensors into Eq. (81c), we determine the coefficients of the linear combination. The 
result is 






AE{d-2) 

\q\ 

327r{d-2)E 



qiqj 



qtqj 



d- 1 



rii 



Aq^ 



-ns 



5'^ + (rf - 3) 



q' 



(82) 



where n^ is defined in Eq. (C3c). Because S^^g and S**-^ vanish, we find that S'4"' = S**^. The 
integral in Eq. (82) is finite and, therefore, we may set d = A. 
The covariant form of the integral 6*4"^ at (i = 4 is then 



S^," 



(47r)2 ~ 



t^u 



1-6' 



m^ 



pf^p'^ 



~5^ 



(83) 
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5. Summary of the NRQCD vertex correction 



Substituting Si - S'^^ in Eqs. (58), (69), (76), and (83) into Eq. (48) and using the 
equations of motion in Eq. (35), we obtain 



A 



NRQCD 



H^ 



NRQCD 



"f^(l+5^) 

47r 


2L{6 


asCp 1 - 5^ 


in 


Att m 


5 



.2 „-^ / 1 TT/i^e-^ 

r r I . log 1 



(I 1 \ 7r2 i-K f I 



3^2 



l + (52 



(84a) 
(84b) 



B. Wave-function Renormalization 



In the Feynman gauge, the self energy of the heavy quark, evaluated at four-momentum 
Pi, is 

P(p,)]nrqcd = -^^-^^f (p + ,,)[(p, + A;)^-^2 + ,,] > (85) 

where m is the mass of the heavy quark and k is the loop momentum, which has been chosen 
to be the momentum of the virtual gluon. In d dimensions, we find that the numerator factor 

reduces to 

2r^ 4r {2 - d){^i + ^) + dm 



P(pi)] 



NRQCD 



-^gtCi 



■h 



(hP' + ie)[{pi + kY — rn? + ie] 
The heavy-quark wave-function renormalization Zq is defined by 



(86) 



[Zq] 



NRQCD 



^_ Pi g[S(pi)]NRQCD 

m dpi 

. Pi '9[S(pi)]nrqcd 



-I -1 



^l=m 



m 



0(a?). 



(87) 



'^i=m 



Differentiating Eq. (86), we find that 



Pi 5P(pi)]nrqcd 



m 



dpi 



-iglCi 



^\=m 



-iglCp 



vAdoDi 



2-d 2[{2 - d) (It + m) + dm]{pi ■ k + m^) 



2-d (2- 



mDoDj 
2m f 1 1 



m 



2m^ 



(88) 



where D^ and D^ are defined in Eq. (34). The expression in Eq. (88) can be written in terms 
of the integrals T02, Tn, T12, T^, Tf^, and T|^, which are defined by 

T. - l-^, (89a) 

(89b) 



rpil. 



fkD'oD'i 
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These integrals are evaluated in Appendix D, and the results are summarized in Eqs. (D7) 
and (D8). The only nonvanishing integral is T^. Hence, 

P^(9[S(pi)]nrqcd 



m 



dp]^ 



AtgiCprn'Tu. 



(90) 



^l=m 



Making use of Eqs. (87), (90), and (D8), we obtain the heavy-quark wave- function renor- 
malization in NRQCD: 



[^q]nrqcd = 1 + ^ f ^ ^ 1 + Oiai). 

2n Veuv em 



(91) 



C. Summary of NRQCD results 



Making use of Eqs. (84) and (91), we find that 
asCp 



G 



NRQCD 



1 



H, 



asCp 1 - 52 



{2[{l+5')L{5)-l](—-^ 
{. ^ eiR eu 



cuv 



NRQCD 



ITT 
1 



4tt m 
Expanding Eq. (28) through order v"^, we obtain 



^: 



QQi 



NRQCD 



r/V^e 



47r 1 3 V^iR euv 






TX^ in f 1 

— + log 

V V V eiR 



TTyU^e 'fe 



2m2 



47r 



— + log 

V V VeiR 



q- 



(92a) 
(92b) 



Sinv 



2m 

V 



(93) 



Comparing Eq. (93) with Eqs. (4.28) and (4.29) of Ref. [9], we find agreement. We have 
also checked Eq. (93) by carrying out a conventional calculation in NRQCD. 
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VII. RESULTS FOR THE SHORT-DISTANCE COEFFICIENTS 

Now we can collect the results of our calculations and obtain the short- distance coeffi- 
cients. By making use of Eqs. (31), (44), and (92), we find that 



^Ga. = ^{2[(i + m.)-i](-i- + log: 



rn? 
+ 6(5"L(5)-4(l + 5")i^(5)-4J>, (94a) 

Aij(^) = ^^ii^^L(5). (94b) 

An m 

As expected, the infrared poles in G*^^) and Gj^j^q^^j-, have canceled in AG*^^-*. Note that 
AG*-^-* and AH^ '' are real and contain only even powers of f = \q\/m. Renormalizing the 
matrix elements in the MS scheme, we obtain 



^C = ^ {2 [(1 + ^')m - 1] log ^ + 66^m - 4(1 + 6^)m - 4} 



(95) 



where now fi is the NRQCD factorization scale. Using Eq. (29), we obtain the short- distance 
coefficients a„ and 6„ : 

of = Sm, (96a) 

do) 3 

bf = -TT-TT- (96d) 



&r = ^^, (96c) 



5 
16m^ 



The results in Eqs. (96a)-(96c) agree with those in Eq. (5.5) of Ref. [6] and those in 
Eqs. (3.13)-(3.20) of Ref. [13]. Using Eqs. (32), (94b), and (95), we obtain the short-distance 
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coefficients 



.(1) 



and 



MS 



^(1) 



J MS 



,(1) 



J MS 



.(1) 



MS 



a« 



J MS 



,(1) 



agCp 1 

47r vn?' 

agCp 1 

47r mf' 

OLsCp 



- + -lo ^ 
9 3 ^"^ 

92 8 

~75 ~ 5 



m^ 



1 -" 



(1) 



&1 



^(1) 



MS 



MS 



MS 



1 / 13744 128, 11^ 
Att m6 V 11025 105 ^ m^ 



a.,C 



s'-^F 



Air m?^ 

— ( n + o log — 



^(1)' 



J MS 



47r m^ \^9 ' 3 

asCp 1/107 9 

47r m6 I 150 5 



m' 



log 






(97a) 
(97b) 

(97c) 
(97d) 
(97e) 
(97f) 

(97g) 



The operators Oao-, Oai, and Obi in Eq- (16) correspond to the operators that were consid- 
ered in Ref. [9], provided that one neglects the gauge fields in the latter operators. There- 



fore, short-distance coefficients [ooIms' ['^i]ms' 
Eq. (4.29) of Ref. [9] as follows: 

[«o]ms 
FiJms 

[^i]ms 



and [b 



iJms 



Cl, 



m^ 



C3, 



2m'^ 



C2- 



are related to the coefficients q in 

(98a) 
(98b) 

(98c) 



Our results for these short-distance coefficients agree with those in Eq. (4.29) of Ref. [9]. 



A. Resummation 

Let us define ratios of the S'-wave QQ operator matrix elements to the S'-wave QQ 
operator matrix element of lowest order in v. 






(99) 



where O^n is defined in Eq. (16a), and we have used the property that the ratios are 
independent of the value of the index i. In Ref. [14], it was shown that these ratios of 
operator matrix elements are related according to a generalized Gremm-Kapustin relation 
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[15]: 

[i<i"nHi^s.)]^=mHes.)Tm- (100) 

This relation holds for the matrix elements in spin-independent-potential models. Hence, 
for each value of n, it holds up to corrections of relative order v'^. 

We can use the relation (100) to resum a class of relativistic corrections to the quarkonium 
electromagnetic current. From Eqs. (20a) and (32), we find that 

oo 

EW+W^]ms)(0|^VI^('5i)) 

n=0 

X (0|OVI-ff('Si)>. (101) 

Because the relation (100) contains corrections of relative order v"^ at each order f ^", the 
resummation in Eq. (101) does not improve the nominal accuracy beyond order f^. The 
resummation might, however, improve the numerical accuracy beyond the accuracy that is 
obtained through order f ^ if the coefficients in the velocity expansion grow rapidly with the 
order in v. In any case, it is interesting to use the resummed result to examine the rate of 
convergence of the velocity expansion. 

B. Numerical results and convergence of the velocity expansion 

Let us evaluate the sums of products of S'-wave short-distance coefficients and operator 
matrix elements, using the relation (100). For {q'^)H{^Si)j we take the central value of the 
J/^ matrix element from Ref. [1]: {q^),j/^ = 0.441 GeV^. Taking rric = 1.5 GeV and setting 
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H = rric, we find that 



Z_^ I ri J MS 
n=0 

t [4"] 

n=0 
n=0 

t [4"] 

n=0 

oo 



MS 



MS 



MS 



(q^)j/v] 
(g^)j/v] 



MS 



MS 



(q^)j/v]5' 



C^sCf 


47r 


CtsCp 


An 


OisCp 


An 


OisCp 


An 


OisCp 



MS 



n=0 



An 



X 7.826, 



X 7.883, 



X 7.872, 



X 7.873. 



(102a) 
(102b) 
(102c) 
(102d) 
(102e) 



In the last hne of Eq. (102), we have used the resummed result in Eq. (101). Taking 
Us = as{2mc) = 0.25, we see that the corrections of order a^f^ and agV^ are 0.5% and 
—0.2%, respectively. These are not very significant at the current level of precision of the 
theory of J/ip decays to a lepton pair. 

As can be seen from Eq. (102), the velocity expansion converges rapidly for approximate 
charmonium matrix elements. In fact, the expressions for AGi^ in Eq. (95) and AH^^' in 
Eq. (94b), taken as functions of f = \q\/m, have finite radii of convergence. The logarithms 
in L{6) [Eq. (41a)] and the Spence functions in K{6) [Eq. (41b)] have branch points at 
6 = ±1, i.e., V = ±oo. The quantity 6 = vj \/\ + v^ has branch points at f = ±i. Therefore, 
the closest singularities to the origin in AG^-^ or AH'^^'' are at f = ±i. Consequently, the 



MS 



radii of convergence of AG^^ and AH^^' as functions off are one. It follows that the velocity 
expansion for the QQ operators is absolutely convergent, provided that the absolute values 
of the operator matrix elements are bounded by a geometric sequence in which the ratio 
between elements of the sequence is less than m'^. 

VIII. CONCLUSIONS 

We have presented a calculation in NRQCD of the order-a^ corrections to the quarkonium 
electromagnetic current. Our calculation gives expressions for the short-distance coefficients 
of all of the QQ NRQCD operators that contain any number of derivatives but no gauge 
fields. Our operators are not gauge invariant, and we evaluate their matrix elements in 
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the Coulomb gauge. Our principal results are given in Eqs. (94b) and (95). The NRQCD 
short- distance coefficients can be obtained, according to Eq. (32), from the Taylor-series 
expansions of AG^ in Eq. (95) and AH^-^^ in Eq. (94b). Our results at relative order v'^ 
agree with those in Ref. [9]. 

Our calculation makes use of a new method for computing, to all orders in v, the one- 
loop NRQCD corrections that enter into the matching of NRQCD to full QCD. In this new 
method, we begin with QCD expressions for the loop integrands. We obtain the NRQCD 
corrections from these QCD expressions by carrying out the integration over the temporal 
component of the loop momentum and then expanding the loop integrands in powers of the 
loop and external momenta divided by the heavy-quark mass m. We carry out this expansion 
before implementing the dimensional regularization. The new approach allows one to avoid 
the daunting task of obtaining NRQCD operators and interactions to all orders in v, along 
with their Born-level short- distance coefficients, and computing their contributions to the 
one-loop corrections. In terms of the total labor involved, the computation of the NRQCD 
corrections to all orders through the new approach is comparable to the calculation of the 
NRQCD corrections at relative order f ^ through conventional NRQCD methods. This new 
method should be applicable to matching calculations for a variety of effective field theories, 
including heavy-quark effective theory and soft-collinear effective theory. 

As we have mentioned, our approach is related to the method of regions [10]. The 
NRQCD corrections in our approach correspond in the method of regions to the sum of 
the contributions from the potential, soft, and ultrasoft regions, i.e., the contribution from 
the small-loop-momentum region [10]. In our approach we have computed the quantities 
AG*-^-* and AH^^^ by subtracting the NRQCD corrections from the full-QCD corrections. 
In the method of regions, AG'-^^ and AH^^^ could, in principle, be computed directly from 
the contribution from the hard region. However, a straightforward computation of the 
contribution from the hard region, carried out by expanding the integrand in powers of 
the small momentum, would yield Taylor-series expansions of AG'-^^ and AH'^'^^ in Eq. (94) 
in powers of 6. It would be nontrivial to sum those expansions to obtain the compact 
expressions in Eq. (94). In contrast, in our approach, expansions of the integrand occur 
only in the NRQCD expressions and lead to very simple series that can be summed at the 
integrand level. Hence, our method may be more efficient than the method of regions for 
computations of short-distance coefficients to all orders in v. Our method is also applicable 
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in the case of a hard-cutoff regulator, such as lattice regularization, while the method of 
regions applies only in the case of dimensional regularization. 

Because we have omitted operators that contain gauge fields, the operators that we 
consider are not the complete set of NRQCD operators that describe the quarkonium elec- 
tromagnetic current. In the Coulomb gauge, the gauge- field operators first enter at relative 
order f^, and so our results cannot be considered to be complete beyond order f ^. However, 
the operators that we consider account for all of the contributions that are contained in 
the Coulomb-gauge wave function of the quarkonium QQ Fock state. The correction to the 
S'-wave component of the electromagnetic current that we find in relative order agV* is only 
about —0.2%, which is not significant at the current level of the precision of the theory of 
J /ip decays to a lepton pair. 

We have examined the convergence of the NRQCD velocity expansion for S'-wave QQ 
operators. In Eq. (102), we give the numerical values for the sums of the first few S'- 
wave contributions to the electromagnetic current and for the sum of all of the S'-wave 
contributions. In these computations, we have made use of the value of the relative-order-f ^ 
J/ip matrix element that is given in Ref. [1] and the approximate relation between operator 
matrix elements in Eq. (100), which holds in spin-independent-potential models [14]. It can 
be seen from Eq. (102) that the velocity expansion converges rapidly in this case. In fact, the 
expressions for AG^ in Eq. (95) and AH^^^ in Eq. (94b), taken as functions of f = \q\/m, 
have radii of convergence one. Therefore, the velocity expansion for the QQ operators is 
absolutely convergent, provided that the absolute values of the operator matrix elements 
are bounded by a geometric sequence in which the ratio between elements of the sequence 
is less than rn^. 
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APPENDIX A: TENSOR-INTEGRAL REDUCTION 



In this Appendix, we describe the tensor-integral reduction that we use to simphfy 
Eq. (36). 

Tensor integrals of rank-1 and -2 that depend on p or on both p and q can be expressed 
in terms of scalar integrals as follows: 

k^f{k,p) = ^ [p-kf{k,p), 
P Jk 



(Ala) 
(Alb) 
k'"f{k,p,q) =p^ I d3{k,p,q)f{k,p,q) + q'' I d4{k,p,q)f{k,p,q), (Ale) 



k^k^f{k,p) = / [d,{k,p)g'^'' + d2{k,p)p^p'']f{k,p), 



k^k''fik,p,q) = g^^ / d,ik,p,q)fik,p,q)+p^p^ / d,ik,p,q)fik,p,q) 

k Jk Jk 



+q''q'' / d7{k,p,q)f{k,p,q) 



+ {py"+P"qn ds{k,p,q)f{k,p,q), 
'k 



(Aid) 



where / is an arbitrary scalar function of the argument four- vectors, li p ■ q = 0, then the 
functions di in Eq. (Al) are given by 



di{k,p 
d2{k,p 
d3{k,p,q 
di{k,p,q 
d^{k,p,q 
dQ{k,p,q 
d7{k,p,q 
d8{k,p,q 



1 

d-l 


e 


1 


{d-l)p^ 


k ■ p 


9 ' 


k ■ q 


q' ' 


1 
d-2 


e 


1 


(d - 2)p2 


1 


(rf-2)g2 


k ■ p k ■ q 



{k-p) 



21 



pZ 
-k^ + d 



2 , Ak.pf 



pz 



{k-pf {k-qf 



pZ gZ 



p. 



p2 



+ {d-l 



q^ 
ik-qf 



q' 



p2 q2 



(A2a) 
(A2b) 
(A2c) 
(A2d) 
(A2e) 
(A2f) 
(A2g) 
(A2h) 
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APPENDIX B: INTEGRALS FOR THE QCD CORRECTIONS 

In this Appendix, we evaluate the integrals in Eq. (38). Throughout this Appendix, we 
neglect expressions of order e or higher. The integrals in Eq. (38) can be expressed in terms 
of elementary integrals /qio, /ho, /on, /-iii, and /m: 

(Bla) 
(Bib) 
(Blc) 
(Bid) 
(Ble) 
Toio, (Blf) 

(Big) 



(B2) 

In deriving Eq. (Bl), we have used the fact that labc = lachi which follows from the sym- 
metry of the integrals under pi ^^ p2 and k — * —k. We have also discarded the scaleless, 
power-divergent integral /loo, which vanishes in dimensional regularization. In deriving the 
expressions for J4 and Jq, we have made a further tensor reduction, using Eq. (Ala), which 
leads to 

/ii-i = — 7 /oio- (B3) 

/qio and /no, which depend only on m^, are given by 

A.o = 7lV(- + '°^^^^ + 2). (B4b) 



Jl-- 


= /oil, 








J2-- 


= /ill. 








h-- 


= /no 


~ /on. 






Ja-- 


1 


2 y''' " 


/-111 - - 

4g2 


^010 \ 


d- 


/ 


J^-- 


= 0, 








Je-- 


= -J4 


+ v^- 




•2m2 


J7-- 


= 0, 








where the scalar integral labc is defined by 










'-'- - I 


1 






.,D-,D\Di 





The scalar integrals /on and /-m can be evaluated by using Feynman parametrization. 
After integrating over k, we obtain 



2 /Anfi^^' '•I 

(4^V^ / .0 



/on = 77^ ( ^ ) r(e) / dz (^2 _ 52 _ ^^ye^ ^35^.) 



/-Ill = tA (^) ^ m (^) P" f dz {z^ -5'- tef-% (B5b) 



(47r)2 V PW V 1 - e / ^0 
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where z = 2x — 1 and x is the original Feynman parameter. Expanding the integrands of 
Eq. (B5) in powers of e, integrating over z, and using Eq. (15), we find that 
i 



-^011 

/-111 



(47r)2 
i m 



1 +log^^^ll^ + 2-26^L{6)+t7T6 



(47r)2 1-6 



euv m 

2 

1-35 



(B6a) 



2A / 1 . log l!i:/f!l^ + 1 ) _ 2^2 ^ 45^^L(5) - 27rt6^ 



^uv 



m^ 



(B6b) 



where L{6) is defined in Eq. (41a) and we have used the following results, which hold for 
0<5 <l: 

f log{z^ -5^- ie) dz = -2 + 25^L{5) + log(l - 5^) - inS, (B7a) 
Jo 
-1 



' z^ - 5"^ - ie) \og{z'^ - 5"^ - ie) dz = - 
3 



+ 46' - A6^L{6) 



+ (1 - 36') log(l - 6') + 2ni6^ 



. (B7b) 



Jill can be evaluated by using Feynman parametrization. After integrating over k, we 

obtain 

•1 



/ill 



4.,^A<r(i + .) r,^^-i-.. 



dz [z' - 6' - ie) 



'l-e 



fB8) 



(47r)2 \ p^ J p^ Jo JO 

where z = 2x—l and the original Feynman parameters are x and y. The infrared divergence 
is isolated in the integral over y: 

1 



dyy 



-l~2e 



(B9) 



'0 2eiR 

The integral over z can be evaluated by expanding the integrand in powers of e. Then, we 
obtain 



/ill 



i 1 - 5M / 1 , infi'e-^ 
— + log 



(47r)2 4m2 I \eiR 



m^ 



2/^(5) + J 



+ 4K{6) 



7i' in , 46' 

-T-J'''^T^6^ 



(BIO) 



where K{6) is defined in Eq. (41b) and we have used the following results, which hold for 
0<6 <1: 

Jq z^ — 6^ — le 26 

' '°«y-f:'^> d. = - iog(i - s^)m - 2im + g + f iog(2«: 

-2"^ — 0^ — «£ 26 6 



(Blla) 
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By making use of Eqs. (Bl), (B4), (B6), and (BIO), we find that 



Ji 



J2 



{Any 



cuv 



+ log 



47r/i^e~^ 



m^ 



2 - 25^L{5) + iTi5 



i 1 - (5M / 1 , 47ru2e-Tfe 
+ log 



(47r)2 4m2 \em 



m^ 



- 2L{6) + 



m 



+ AK{6) 



TX 



m , 45^ 
log- 



b 6 ^1-52 



J. 



J4 



{Any 



{Any 



25^L{5) - iTiS 



+ Ji 



Jps = 0, 



Jj = 0. 



The results for J1-J4 in Eq. (B12) agree with those in Ref. [16]. 



APPENDIX C: INTEGRALS FOR THE NRQCD CORRECTIONS 



(B12a) 



(B12b) 

(B12c) 

(B12d) 
(B12e) 
(B12f) 
(B12g) 



Here, we tabulate some integrals that are useful in computing the NRQCD corrections. 
In dimensional regularization, scaleless, power-divergent integrals vanish: 

1 



k ll^l 







for n 7^ 3. The only scaleless logarithmically divergent integral that we encounter is 

[ — = -(—--)■ 



no 



(CI) 



(C2) 



There are a few integrals that depend on q that appear in the evaluations of the Si in 
Eq. (49): 



111 



^2 



1 



ir\l\ 



k"^ + 2k ■ q -ie An 
1 



(C3a) 



fc fc2(fc2 ^2k ■ q-ie 
i ( 1 



nz 



167r|qr| ^eiR 

fc2 



log 



TT/i^e t: 



in 



k"^ + 2k- q-ie 2n 



i^m ■ 



(C3b) 
(C3c) 
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We also make use of the angular averages: 



im 1 ,„^ (E^ I , ,(,.)_ ,c4a, 



kE±q-k 2|g| \E - \q\^ j^ 



fik') T—^ + 



kE^-{q-kf 2E Jk \E + q-k E-qk 

1 .^^jE+\q\\ /•^,,2 



log ^r^ /(fc^), (C4b) 



2E\q\ ° \E 
where f{k'^) is any function of k"^. 



APPENDIX D: EVALUATION OF THE INTEGRALS FOR [Zq]nrqcd 

In this Appendix, we evaluate the integrals T02, Tn, T12, T^, Tf^, and T^^, which enter 
into the calculation of [2'q]nrqcd and are defined in Eq. (89). We make use of the same 
strategy that we used in evaluating the Si integrals in Sec. VI, except that we carry out the 
evaluation in the rest frame of the heavy quark, pi = {m, 0), where the expressions become 
compact. The change of frame shifts momenta by an amount of order mv. Therefore, the 
NRQCD expansion in powers of the external momentum divided by m remains valid. In the 
heavy-quark rest frame, the gluon- and quark-propagator denominators are 

[^0] rest = (^° + 1^1 - ^^)(^° - \k\ + ie), (Dla) 



[Di] ^^^^ = ( A;° + m + Vm'^ + k^ - ie){k^ + m - Vm'^ + k^ + ie). (Dlb) 

We evaluate the k^ integrals by using contour integration, closing the contour in the upper 
half-plane in every case. We denote the contributions of gluon and quark poles by subscripts 
g and Q, respectively. 
The integral T02Q is 



The integral Tu yields 



^ 02Q ^ 






^^'^—irnti^['-v^;^TW^- '°''=' 
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The integral T12 yields 



Ti 



129 



6m^ 



\k\ 



im" 



■no, 



Ti 



12Q 



6m- 



.4r 1 

Jk {ni? + 



fc2)3/2' 



(D4a) 
(D4b) 



where Uq is defined in Eq. (C2). 

In the cases of the integrals T^, T^^, and T{^, the integrand of the temporal component 
T°j is identical to that of Tab-, except that the integrand in T^^^ contains an additional factor 
of A;°. Integrating over kP , we obtain 



rpG 

^02Q 



im A 1 



fc2)3/2' 






i<y 



Ik 



rpO 
-^1 



IQ 



^123 



rriO 

^12Q 



^i- 



1 



i f 1 



5m^ 



\k\ 



6m- 



4\ 

Jk k 



1- 



m 



m^ 



fc2)3/2 



(D5a) 
(D5b) 
(D5c) 
(D5d) 
(D5e) 



For the spatial component T^^,, the integrand is identical to the integrand in Tab, except that 
the integrand in T^^ contains an additional factor k^. By making use of Eqs. (D2)-(D4), we 
find that 

4 7fc(^' + fc')^/'' 



T 



02Q 



T^ = / 

'^' Am A fc2 ' 



k' 



Ti 



m 



IQ 



T 



4m f^k^ \ Vm2 + k^ 
i f k 



I2g 



6m? ./fc \k\'^' 



T 



12Q 



5m' 



Jk 



t 



m^ 



^2)3/2- 



(D6a) 
(D6b) 
(D6c) 
(D6d) 
(D6e) 



Expanding the integrands in Eqs. (D2)-(D6) in powers of k"^ /ir?, we find that all of the 
terms in the expansions yield scaleless, power-divergent integrals, with the exception of the 
integral T123 in Eq. (D4a). Therefore, 



T 



11 



^02 



-02 



"'M T^M n 



-^11 ~ -^] 



12 



(D7) 
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and 

Ti2 = - ^^ : . {- ^ 1 , (D8) 



327r2m2 ^euv em 
where we have used Eq. (C2). 
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